Probing correlated phases of bosons in optical lattices via trap squeezing 
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We theoretically analyze the response properties of ultracold bosons in optical lattices to the 
static variation of the trapping potential. We show that, upon an increase of such potential (trap 
squeezing), the density variations in a central region, with linear size of > 10 wavelengths, reflect 
that of the bulk system upon changing the chemical potential: hence measuring the density vari- 
ations gives direct access to the bulk compressibility. When combined with standard time-of-flight 
measurements, this approach has the potential of unambiguously detecting the appearence of the 
most fundamental phases realized by bosons in optical lattices, with or without further external 
potentials: superfluid, Mott insulator, band insulator and Bose glass. 

PACS numbers: 03.75.Lm, 71. 23. Ft, 68.65.Cd, 72.15.Rn 



Ultracold gases in optical lattices offer the unique op- 
portunity of literally implementing fundamental lattice 
models of strongly correlated quantum many-body sys- 
tems, either bosonic or fermionic, traditionally consid- 
ered as "toy" models for the description of complex con- 
densed matter systems 0, 0| ■ In the particular case of ul- 
tracold bosons realizing the Bose-Hubbard (BH) model, 
recent experimental developments have led to the spec- 
tacular demonstration of the Mott insulating (MI) phase 
with controllable filling [H, 0, 0| , and even more recent 
developments in laser trapping offers the possibility of 
realizing further fundamental insulating phases, such as 
a band insulator (BI) in a commensurate superlattice Q 
or a Bose glass (BG) in an incommensurate superlattice 
or in a laser-speckle potential @, H| ■ 

Two main technical aspects limit to date the possibility 
of the experiments to retrieve full information on the true 
bulk behavior of the model Hamiltonian implemented in 
the system. One aspect is the presence of a parabolic 
trapping potential which imposes a spatial variation of 
the filling and hence a spatial modulation of the local 
behavior exhibited by the system. A second aspect is 
represented by the typical measurements performed on 
the system. The detection of strongly correlated phases 
is typically based on the measurement of correlation func- 
tions (phase correlations via time-of-flight measurements 
3| and density-correlations via noise-correlation analysis 
9]) and on site-occupation statistics 0, [l(J. Lattice- 
modulation spectroscopy offers the possibility of mea- 
suring the dynamic structure factor at zero transferred 
momentum [8|, [ll| but it has the drawback of probing 



the global response of the inhomogeneous system, and 
of being subject to a low-energy cutoff imposed by the 
duration of the experiment. This aspect prevents e.g. 
the unambiguous observation of the BG, which does not 
have a special signature in correlation functions, but it 
is unambiguously marked by the absence of a gap in the 
excitation spectrum. 

The purpose of this paper is to propose a technique 
- trap- squeezing spectroscopy - which circumvents these 



two limitations at once, taking advantage of the parabolic 
trapping to extract the bulk behavior of the model im- 
plemented in the system and in particular its lowest 
particle-hole excitation energy. Two fundamental obser- 
vations are at the basis of this proposal. On the one 
hand, the slowly varying nature of the parabolic poten- 
tial guarantees the validity of the local-density approx- 
imation (LDA) [ill. HH, particularly close to the poten- 
tial minimum. Hence the density around the trap center 
mimics the behavior of the bulk system over an extended 
region of space of linear size of several (>10) lattice spac- 
ings, and consequently the average density in this region 
can be accessed via laser microscopy [13]. On the other 
hand, this central average density can be controlled via 
the trapping potential in very much the same way as 
the chemical potential controls the density of a bulk sys- 
tem in the grand-canonical ensemble. In particular the 
trapping potential is by far the lowest-energy potential 
to which the system is coupled (with trapping frequen- 
cies as low as ~10 Hz), and measuring the response of 
the central density to small variations of such potential 
allows to directly probe the low-energy response of the 
bulk system. 

We theoretically investigate the trap-squeezing spec- 
troscopy in the one-dimensional (Id) BH model in a 
parabolic potential plus an external superlattice poten- 
tial, n(j,u,v 2 ,v t ) = n (j,u,v 2 ) + v t J2S 

where 

Ti-o = ^2 -J (b\b l+1 + h.c.) + — n,(ni - 1) + V 2 g t n 

i 

(1) 

Here <?, = gi(a, 4>) — cos 2 (27ra i + <fi) ~ 1/2 is a one-color 
superlattice potential. J and U are experimentally con- 
trolled via the height of the primary optical lattice, while 
V2 is controlled via the height of a secondary optical lat- 
tice [1, Q. In the following we make the fundamental 
assumption that Vt can be varied independently of the 
other parameters, which is possible by applying an ex- 
tra dipolar trap to the system created by an additional 
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running laser wave. 

We study the above model via Stochastic Series Ex- 
pansion quantum Monte Carlo at low temperatures 
(capturing the T = behavior) and in the grand- 
canonical ensemble, namely we simulate the Hamilto- 
nian TL^ = TL — \i J^. m where the chemical potential 
[i is fine tuned to get the desired average number of par- 
ticles (N) , and in this way it becomes a function of the 
other Hamiltonian parameters fi = fi(Vt, N, J, U, V2). Ac- 
cording to LDA, the average density at the center of the 
trap ric —■ 1/|C| J2 ieC {ni) (where the region C will be 
defined later) reproduces closely that of a homogeneous 
system (V* = 0) at a chemical potential /1. Hence con- 
trolling /i via one of the other parameters Vt, N, J, U, 
and V2, allows to control nc- In particular, if fj, is con- 
trolled by changing Vt, namely by trap squeezing, while 
holding all the other parameters fixed, one has access to 
the compressibility for the bulk Hamiltonian Ho(J, U, V%), 
estimated via k — dnc/d^i. 
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FIG. 1: Scaling of the chemical potential n for the Id Bose- 
Hubbard model in a parabolic trap (a) , and in a trap plus an 
incommensurate superlattice potential with strength V2 = U 
(b). fi,Vt,V2 and U are here reported in units of J. The 
dashed lines correspond to a linear fit to the lowest-(t// J) 
data, rescaled by a factor (U/ J) 2//3 to compare with the other 
data sets. The inset shows the collapse of all curves over a 
universal curve / M (see text). 

The control on the chemical potential \x via trap 
squeezing requires the detailed knowledge of the func- 
tion /j, = /J,(Vt, N, J, U, V2). Such a function can be ac- 
curately sampled via quantum Monte Carlo, given that 
its values are the result of the fine-tuning procedure of 
the chemical potential required to achieve a desired av- 
erage N, Fig. Q] shows [i = n(Vt , N, ...) for different cases 



of the BH model without external potentials, V2 = 0, 
and for an applied incommensurate superlattice potential 
with strength V2 = U and incommensurability parameter 
a = 0.7714... identical to that of the experiment of Ref.d. 
The behavior of /i for a commensurate superlattice with 
a = 3/4 is found to be nearly identical to that of the 
incommensurate case. In absence of a superlattice, and 
for weakly interacting bosons, Thomas-Fermi (TF) the- 
ory [15| would predict the following scaling for the chemi- 
cal potential in (i-dimensions: fi ~ (N [/) 2 /( 2+d ' yd/( 2 + d ) 
which gives fi ~ (N 10 2/3 V t 1/3 for d= 1. Fig. Q] shows 
that, at fixed U and for all the cases considered, /z is a ho- 
mogeneous function of the combination x = N 2 / 3 V t ^ 3 ; 
in particular, even for large U /J ratios it suprisingly ver- 
ifies the TF prediction of a linear dependence on x; sig- 
nificant deviations are observed only in the low-density 
and hi gh-( ?7/J) case, where the hardcore boson regime 
sets in [16| . Moreover, for sufficiently low U/J the data 
for different f/'s collapse on the same universal curve 

t ,= (U/jf/ 3 U(x) (2) 

where / M is essentially a straight line. A linear fit for the 
lowest- (£/ 1 J) data gives = -1.225(13) + 0.817(2)x 

for all the cases considered, namely in presence or in 
absence of a superlattice. We notice that the crude 
TF prediction would give f^x) = [(d + 2)r(eJ/2 + 
l)/(2ir d / 2 )] 2 /( 2+d ' > x, which is well off the numerical data. 
For higher U/J we observe that the (U/J) 2 / 3 scaling is 
still essentially obeyed by the derivative d\x/dx but not 
by the intercept fi(x = 0); hence the slight disagreement 
between the scaling curve and the data for U / J = 20 in 
Fig. [Ha) and U/J = 30 in Fig. \V(h). 

Therefore we obtain a universal prediction for the de- 
pendence of the effective chemical potential in the center 
of the trap on the experimentally controllable parame- 
ters J, U, V t and N for a large range of their values, 
and for the extreme case of d — I where the applicabil- 
ity of mean- field theory is in doubt. Similar results are 
then expected to hold a fortiori for the cases d — 2,3. 
Hence we can firmly conclude that the chemical poten- 
tial in the center of the trap represents a well controlled 
experimental parameter. 

Armed with this prediction, we can then move on to 
simulate the outcome of a trap-squeezing experiment, 
where the central density nc is monitored as a function of 
the trapping potential Vt. We start from the case of the 
Id BH model without any superlattice, for which we con- 
sider a boson number (N) — 100 in a variable-frequency 
trap and with fixed repulsion U/J = 20. Fig. [2] shows 
the evolution of the central density nc averaged over a 
region C containing 10 — 20 sites as a function of the 
chemical potential fJ>(Vt), and compared with the data 
for the bulk system. It is evident that, for a sufficiently 
low \i (namely for sufficiently low Vt), the bulk density 
curve is very well reproduced (in this case for n c < 2). 
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FIG. 2: Central density nc and global coherent fraction Uk=o 
for the Id Bose-Hubbard model in a trap as a function of 
the effective chemical potential required to mantain N = 100 
bosons in the system. Here U/J = 20. The 'bulk system' 
data are obtained via a simulation on a homogeneous system 
with negligible finite-size effects and in the grand-canonical 
ensemble with chemical potential /i. 

The deviation of nc from the bulk value reveals that the 
truly homogeneous region in the trap center has become 
smaller than the C region, a fact that can be simply 
cured by increasing the number of particles and decreas- 
ing the trapping potential so as to leave /i ~ V^ 3 N 2 / 3 
fixed. The succession of incompressible plateau regions 
at integer filling and compressible regions in the n c (p) 
curve marks the alternation between incoherent MI and 
coherent superfluid (SF) behavior, as also revealed by 
the (global) coherent fraction nk=o — (l/-^0 Xh? W&j)- 
Remarkably, when the effective chemical potential in the 
trap center overcomes the Mott gap, a few particles can 
be transfered from the wings to the center into a locally 
SF state, and this gives rise to a violent increase in the 
coherent fraction with a very sharp kink. The width of 
the integer-filling plateaus corresponds to that of the MI 
lobes in the phase diagram of the Id BH model: hence 
this kind of measurement allows to reconstruct that phase 
diagram with high accuracy, and to extract the particle 
(hole) gap at any point as the minimal chemical poten- 
tial variation required to increase (decrease) the density. 
In particular, trap squeezing probes the density-driven 
transition from MI to SF, which is in a different univer- 
sality class 1] with respect to the transition driven by the 
J/U ratio and probed so far in experiments [H, EH H3|- 
Moreover we emphasize the high tolerance of the method 
to the variation of the size of C, which corresponds to 
the size of the focus of the imaging laser. 

Having shown that trap squeezing allows to reconstruct 
the phase diagram of the bulk Bose-Hubbard model, we 
generalize this approach to probe other phases of corre- 
lated bosons in an optical lattice. To this end we con- 
sider N = 100 trapped bosons in an additional com- 
mensurate superlattice potential (l7| with a = 3/4, fixed 
phase 4> = 0, and strength V-x = U = 20J, such that 
it overcomes the MI gap and hence it removes the MI 
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FIG. 3: Central density and global coherent fraction for the 
Id Bose-Hubbard model in a trap and in a commensurate 
superlattice (V 2 = U = 20 J, a = 3/4, <f> = 0). All symbols 
and notation as in Fig. [2] Notice that the deviation of the 
data for C = 10 sites from the bulk ones is due to the fact that 
the C region does not contain an integer number of periods 
of the superlattice potential. 

phase; the insulating phase which is left for large U / J 
is a BI with fractional, commensurate fillings (2n + l)/4 
(n — 0, 1, ...). Fig. [3] shows the alternation of phases in 
the center of the trap under trap squeezing as revealed 
by the central density, and compared to the bulk result; 
similarly to the MI-SF transition, the BI-SF alternation 
is clearly evidenced. The density plateaus correspond 
to the formation of incompressible BI region in the trap 
center, an event associated with a significant lowering of 
the global coherence in the system, as shown by the nk=o 
curve; the coherence is suddenly increased when the BI 
gap is overcome by the chemical potential and particles 
are transfered into a locally SF state in the center. 
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FIG. 4: Central density and global coherent fraction for the 
Id Bose-Hubbard model in a trap and in an incommensurate 
superlattice (V2 = U = 20J, a — 0.7714..); (.■.)$ denotes the 
average over fluctuations of the spatial phase <j>. The boxes 
mark some relevant extended regions exhibiting Bose glass 
(BG) behavior. All other symbols and notations as in Fig. [2] 

The situation changes drastically when tuning slightly 
the superlattice parameter from the commensurate value 
a = 3/4 to the incommensurate value a = 0.7714.. real- 
ized in recent experiments In this case, for a strong 
superlattice V2 = U and for small J/U the ground state 
of the system changes from SF to incompressible incom- 
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mensurate band insulator (IBI) and to compressible BG 
upon changing the chemical potential [13| . Fig. [4] shows 
the variation under trap squeezing for the central density 
averaged over random fluctuations of the spatial phase, 
( n c)(f>- This average is intrinsic in current experimen- 
tal setups, where the phase <f) can change from shot to 
shot, and it is essential for the central region of the trap 
to sample the full statistics of the quasi-periodic poten- 
tial and hence to mimic the bulk behavior of the system 
[l3j ] . Indeed we observe that (nc)<p reproduces very well 
the bulk behavior for low enough density. In striking con- 
trast to the previous two cases of no superlattice and of a 
commensurate superlattice, the {nc)rf, curve exhibits ex- 
tended compressible regions for which the coherent frac- 
tion does not vary upon changing the chemical potential. 
This corresponds to transfer of particles at no energy cost 
from the wings to the center of the trap into localized 
states which do not contribute to the coherent fraction 
of the system: this fact provides smoking-gun evidence 
for the appearence of a BG state in the center of the trap 
[lij ]. Moreover the joint information coming from the 
central density and the global coherent fraction enables 
to experimentally probe the incompressible IBI behavior 
and the compressible SF behavior. 

In summary, we have proposed an experimental 
method (trap squeezing spectroscopy) to directly extract 
bulk properties of strongly correlated bosons from mea- 
surements on a trapped system - a fundamental require- 
ment in the future perspective of quantum simulations of 
complex quantum systems realized with cold atoms. The 
method relies on a simple, universal relationship between 
the trapping potential and the effective chemical poten- 
tial for the particles in the trap center, which we numer- 
ically elucidate in the case of the Bose-Hubbard model 
realized in optical lattices. Measuring the response of the 
central density in the trap to the variation of the trap- 
ping potential provides direct access to the compressibil- 
ity of the infinite system, a piece of information which is 
not directly accessible to current experimental setups and 
which is crucial to extract the energy gap over the ground 
state of the Hamiltonian implemented in the system. The 
method offers the possibility to extract the phase diagram 
of the Bose-Hubbard model with high resolution. Most 
remarkably, the joint measurement of the compressibility 
and of the coherent fraction (obtained via time-of-flight 
techniques |3|) provides clear evidence for the realization 
of a Bose-glass state in the center of the trap. We have 
demonstrated this property in the case of an incommen- 
surate superlattice as recently realized in experiments Q , 
although the same technique can be applied to different 
realizations of random or pseudo-random potentials ■ 

From the experimental point of view, this method re- 
quires the application of an extra dipolar trap whose 
strength can be controlled independently of that of the 
optical lattice, and the measurement of the optical depth 
of the cloud over a region of order ~ 10 wavelengths 



(~ 10/im) of the optical lattice in all three spatial di- 
rections. The measurements of the central density and 
of the coherent fraction cannot typically be performed in 
the same shot, so that special care is needed in mantain- 
ing the number of particles N fixed from shot to shot 
to achieve the same experimental conditions. This can 
be typically achieved by post-selecting only those mea- 
surements with the same total ./V in the trap. Useful 
discussions with N. Bar-Gill, L. Fallani, C. Fort, and M. 
Rigol are gratefully acknowledged. 
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